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We  consider  the  Euler  equations  for  a  perfect  fluid  in  a  flat- 
bottomed  canal  in  the  time-dependent  case.  A  formal  expansion  procedure 
for  small  amplitude,  long  waves  analogous  to  that  of  Friedrichs  and  Hyers 
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tion  (KdV  for  short)  for  the  lowest  order  term.  The  higher  order  terms  in 
the  expansion  satisfy  the  inhomogeneous  version  of  the  linearized  KdV 
equation. 

Of  particular  interest  to  us  are  those  solutions  of  the  KdV  equation 
called  N-solitons,  which  asymptotically  separate  into  N  travelling  waves 
with  distinct  speeds.  Using  certain  facts  about  the  linearized  KdV 
equation  and  some  properties  of  the  N-solitons,  we  prove  that  the  next 
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conditions  and  a  symmetry  requirement.  This  solution  behaves  like  an 
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SIGNIFICANCE  AND  EXPLANATION 


T 


The  Euler  equations  for  fluid  flow  in  a  flat-bottomed  canal  may  be 
approximated,  in  the  long  wavelength,  small  amplitude  limit,  by  the 
Korteweg-de  Vries  equation  (KdV  for  short) .  We  study  the  nature  of  this 
approximation  and  attempt  to  justify  it  mathematically.  The  Friedrichs 
and  Hyers  proof  [9]  of  the  existence  of  the  solitary  wave  provides  such 
a  justification  for  travelling  wave  solutions  of  the  Euler  equations. 

In  fact,  they  gave  a  formal  expansion  procedure  for  the  solution  which 
yields  the  (time-independent)  KdV  equation  for  the  leading  term. 

In  this  paper,  time  dependent  solutions  of  the  Euler  equations 
are  examined  in  the  same  long  wavelength  limit.  A  formal  expansion 
procedure  analogous  to  that  of  Friedrichs  and  Hyers  is  developed,  and 
the  fulL  KdV  equation  is  derived  for  the  leading  term.  Choosing  the 
so-called  N-soliton  solution  of  the  KdV  equation,  we  seek  solutions  to 
the  Euler  equations  with  similar  properties.  The  N-solitons  are  non¬ 
linear  superpositions  of  N  solitary  waves  with  distinct  speeds,  and 
thus,  for  large  positive  or  negative  times,  they  decouple  into  N 
travelling  waves.  By  analyzing  the  linearized  KdV  equation,  we  show 


that  with  suitable  boundary  conditions  and  a  symmetry  condition,  the  next 

term  in  the  formal  expansion  is  unique  and  resembles  an  N  soliton.  The  / 

/  ^  g 

effect  of  the  first-order  correction  to  the  KdV  equation  is  merely  to  /  ~j ,  L>  ■>, 

alter  the  shape  of  the  waves  slightly;  their  speeds  and  phases  remain,,  v  A- 

/  ■  '  t>_  '•'o’ 


the  same. 
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A  JUSTIFICATION  OF  THE  KdV  APPROXIMATION 

TO  FIRST  ORDER  IN  THE  CASE  OF  N-SOLITCN  WATER  WAVES  IN  A  CANAL 

Robert  L.  Sachs 

1.  INTRODUCTION 

The  Korteweg-de  Vries  equation  (KdV  for  short)  was  originally  derived 
in  1895  as  an  approximation  for  fluid  flow  in  a  flat-bottomed  canal  [14] . 
This  non-linear  evolution  equation  for  a  function  of  one  space  variable 
has  the  rather  remarkable  property,  discovered  by  Gardner,  Greene,  Kruskal, 
and  Miura  [10] ,  that  it  may  be  solved  more  or  less  exactly.  In  fact,  a 
Hamiltonian  structure  can  be  introduced  and  the  KdV  equation  may  be 
regarded  as  a  completely  integrable  Hamiltonian  system.  One  very 
interesting  class  of  solutions  is  the  set  of  so-called  N-solitons. 

These  solutions  behave,  for  large  positive  and  negative  times,  like  N 
exponentially  decreasing  'bumps'  moving  at  distinct  speeds.  A  natural 
question  to  ask  is  whether  such  "N-tuple  waves"  exist  for  the  full  set  of 
Euler  equations  governing  the  fluid  flow  in  a  canal. 

For  N  =  1,  such  wave  solutions,  known  as  solitary  waves,  do  in 
fact  exist  [3,4,9].  In  [9],  Friedrichs  and  Hyers  qave  a  formal  expansion 
procedure  for  the  Euler  equations  in  which  a  time-independent  form  of  the 
KdV  equation  arose  as  the  equation  satisfied  by  the  lowest  order  term. 

The  higher  order  terms  of  their  expansion  satisfied  the  inhomogeneous 
form  of  the  linearization  of  the  non-linear  ordinary  differential  equation 
for  the  leading  term.  With  a  symmetry  condition  added  to  the  requirt nor.t 
of  exponential  decay, this  equation  could  be  solved  uniquely.  After  re¬ 
formulating  the  problem,  the  convergence  of  this  formal  solution  wu; 
shown  by  the  implicit  function  theorem.  Later  Beale  [4|  impl  i  1  iel  the 
argument  by  using  a  generalized  implicit  function  theorem  due  to  :'ehnd“! 
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[251 .  In  both  of  these  approaches,  the  time- independent  nature  of  the 
problem  is  relied  upon  from  the  beginning. 

In  attempting  to  generalize  these  results  to  N-solitons  for  N  S  2, 
the  problem  becomes  unavoidably  time-dependent.  An  essentially  trivial 
step  in  both  approaches  to  the  solitary  wave  problem,  namely  inverting  the 
linearized  KdV  operator,  now  becomes  a  serious  difficulty.  Constructing 
a  formal  solution  which  behaves  like  a  N-soliton  requires  solving  the  in¬ 
homogeneous  linearized  KdV  equation  with  prescribed  asymptotic  behavior. 

We  do  this  for  the  first  order  correction  term  by  using  the  explicit  form 
of  the  inhomogeneous  term.  For  higher  order  corrections,  the  existence 
of  some  solution  is  guaranteed  by  Duhamel's  principle  and  the  solvability 
of  the  Cauchy  problem  for  the  linearized  KdV  equation  [19] .  However,  in 
such  an  approach,  initial  values  (say  at  t  =  0)  ’parameterize'  the  set 
of  all  solutions  and  we  cannot  as  yet  single  out  those  solutions  with  the 
desired  asymptotic  behavior. 

In  this  paper  we  present  the  following  results: 

(i)  The  time-dependent  analogue  of  the  formal  expansion  of  Friedrichs- 
Hyers  [91  is  developed.  For  perturbations  of  a  steady  horizontal 
flow  with  Froude  number  near  1  which  are  of  small  amplitude, 
long  wavelength,  and  slow  time  variation,  we  consider  a  formal 
power  series  solution  of  the  Euler  equations.  The  small  per- 
ameter  e  is  related  to  the  Froude  number.  As  in  [9] ,  the 
lending  term  satisfies  the  KdV  equation  and  the  higher  order 
terms  satisfy  the  inhomogeneous  linearized  Kc*V  equation.  How¬ 
ever,  in  this  case,  both  of  these  equations  are  time-dependent. 

(i i )  Using  results  on  the  solvability  of  the  Cauchy  problem  for  the 

linearized  KdV  equation  [19!  and  certain  facts  about  N-solitons, 
we  analyze  the  first  order  term  completely.  In  particular,  we 
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show  that  tnis  term  is  uniquely  determined  by  the  following 
conditions : 

(a)  (symmetry)  u(x,t)  =  u(-x,-t) 

(b)  (asymptotic  decay  in  moving  frames) 

u(ct  +  £,t)  -*  0  exponentially  fast  as 

t  +»  for  f  fixed  unless  c  =  c 

3 

j  =  1,...N  where  {c.}  are  the  N 

3 

soliton  speeds 

(c)  (asymptotic  shape) 

lim  u(c.t  +  £  ,t)  is  an  exponentially 

t-Xo  3 

decreasing  function  of  £. 

This  is  the  sense  in  which  we  use  the  term  justification  in  the  title  of 
this  paper.  The  first  order  correction  to  the  KdV  N-soliton,  as  chosen 
above,  does  not  alter  any  of  the  essential  features  of  the  solution. 
After  a  long  time,  the  water  wave  decomposes  into  N  travelling  waves 
with  distinct  speeds,  each  of  which  is  exponentially  decreasing  in 
space  when  viewed  from  the  appropriate  moving  frame  of  reference. 

Section  2  contains  the  time-dependent  analogue  of  the  formal 
expansion  of  Friedrichs  and  Hyers  [9]  as  well  as  the  mapping  formulation 
of  Beale  [4] .  For  the  latter  set-up,  invertibility  of  the  linearized 
mapping  at  c  =  0  is  shown  in  the  cormal  sense  provided  the  linearized 
KdV  operator  is  invertible.  The  basic  facts  concerning  the  Cauchy 
problem  for  the  linearized  KdV  equation  are  presented  in  Section  3. 
Explicit  solvability  for  this  problem  is  related  to  the  so-called 
inverse  scattering  method  for  f Giving  the  KdV  equation  [10,  19].  Using 
certain  facts  about  N-solitons,  which  we  present  in  the  Appendix,  and 
the  particular  terms  arising  in  the  expansion  of  Section  2,  the  first- 
order  correction  to  the  N-soliton  is  analyzed  in  Section  4. 
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2.  THE  EULER  EQUATIONS,  THE  KdV  LIMIT,  AND  A  FORMAL  EXPANSION 


In  dimensionless  variables,  the  Euler  equations  for  a  perfect 
fluid  in  a  two-dimensional,  flat-bottomed  domain  D,  with  a  free  boundary 
y  =  r  (t,x)  as  upper  surface,  subject  only  to  gravitational  acceleration 
g  are  (cf.  Stoker  [21]): 


(?.l)  < 


(i)  +  <t>  =  0  in  D  =  {  (x,y)  :  0  <  y  <  T(t,x)} 

xx  yy 


(ii)  4>y  =  0 


along  y  =  0 


12  2 

(iii)  +  —  (4>x  +  4>y)  +  yy  =  constant 

along  y  =  r (t,x) 


(iv)  F.  +  <j>  •  T  —  4>  =  0  along  y  =  F  (t,x) 

t  x  x  y 


-  gh 

where  <J>  =  <})(x,y,t)  is  the  velocity  potential  and  y  =  where  h  is 

IT 

the  length  scaling  and  U  is  the  velocity  scaling. 

-1/2 

V  is  called  the  Froude  number  or  reduced  depth  and  is  a  parameter  of 

the  problem.  The  linear  theory  of  water  waves  [21]  predicts  y  =  1, 
while  the  existence  of  solitary  waves  occurs  for  y  <  1  but  sufficiently 
close  to  1.  From  now  on,  we  assume 

(2.2)  0  <  1  -  y<<y<l  and  in  fact,  we  define  a  small  parameter  e 

,  ,  -3e 

by  the  equation:  y  =  e 

In  this  section,  we  will  consider  flows  which  are  very  nearly  the 
trivial  flow  of  constant  horizontal  speed  1  given  by  the  solution 
$  -  x,  F  =  1,  y  =  1  of  (2.1).  Introducing  auxiliary  variables  £,',n’ 
which  vary  over  a  fixed  horizontal  strip  0  <  nr  <  1,  we  may  eliminate 
the  unknown  free  surface  at  the  expense  of  defining  x,y  as  functions 
of  f, ' ,  n',  t.  In  steady  flow  problems,  t  does  not  appear  and  C*  +  in' 
is  usually  the  complex  potential  function,  but  for  time-dependent 
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I 


t 


problems,  we  express  both  the  potential  function  <p  and  the  physical 

coordinates  x,y  in  terms  of  £ ' ,  n*  and  t.  Provided  the  mapping 

(f,,h')|->p  (x,y)  is  invertible  for  every  t,  solving  the  problem  in  the 

'  ,  n*  plane  is  equivalent  to  solving  the  original  system  in  the  x,y 

plane.  In  the  neighborhood  of  the  trivial  horizontal  flow,  this  mapping 

is  roughly  the  identity  map,  hence  it  will  be  invertible. 

After  expressing  the  problem  in  these  new  independent  variables,  a 

new  dependent  complex  variable,  A’  -  iO * ,  defined  as  the  logarithm  of 

the  complex  velocity  W  (W  :  $  -  i<J>  )  ,  will  be  introduced.  By  different- 

x  y 

iating  with  respect  to  f '  along  n*  =  1,  ^  is  eliminated  and  a  new 

system  of  equations  for  x,  y,  A',  O'  is  obtained.  Defining  a  small 
1/2 

parameter  a  =  t  ,  we  rescale  the  independent  variables  £  '  ,  i',  t  and 

the  small  dependent  variables  A’,  O',  x'  x-f, ',  y'  B  y-n '  .  The  system  (2.1) 

in  the  rescaled  variables,  £,  r,  ;  and  x,  y,  A,  9  respectively,  becomes: 


“  J>  — 


Thu:'  ’  derivatives  of  f  are  expressible  in  terms  of  w  and  z. 
Substitution  in  (2.1)  and  differentiation  with  respect  to  £*  along  n'  =  1 
yields  a  system  with  w,z  as  dependent  variables,  namely: 


-6- 


(2.5) 


(i) 

(ii) 

(iii ) 

(iv) 


w(C',t)  and  z(?’,t)  are  holomorphic  functions  of 
in  0  <  Im  ?'  <  1 

Im  w  =  0  ;  Im  z  =  0  along  n'  =  0 

X  2 

Re  (w^,  -  w^,zfc)  +  —  ( J w |  )^f  +  y  Im  (z^f)  =0 
along  n ’  =  1 

Im  (z  /z  ,)  +  Im  (w/z. , )  =  0  along  n 1  =1 

(This  last  condition  comes  from  the  relations  T  =  y  ,/x  , 

X  s 

rt  =  yt  -  vw  on  n'  =  1*) 


It  is  convenient  to  replace  w  by  X'  -  i0*  ,  defined  by  the  relation 


(2.6) 


X»-i0* 

w  =  e 


This  substitution  was  introduced  by  Levi-Civita  [16]  in  the  periodic  case 
of  infinite  depth;  it  has  the  virtues  of  simplifying  the  |w|  differentiation, 
ensuring  w  J-  0  for  any  solution,  and  making  X '  -  i0 f  =  0  the  trivial 
flow.  Upon  substitution,  we  obtain: 

(i)  z(c',t),  (X*-i0 ' ) (£' ,t)  are  holomorphic  in 

for  0  <  Im  <  1 

(ii)  0*  =  0,  y  =  0  along  n'  =  0 

(2.7)  •  , t 

(iii)  Re  (e  ~  [A£  -  i0£)z  f  -  (A£f  -  i0*,)ztJ) 

2X  f 

+  e  X’,  +  Y  Im  z^,  =  0  along  n*  =  1 

X 1  — i0  1  — 

(iv)  Im  (zt/z^,)  +  Im  (e  /z^ f )  =  0  along  n'  =  1 
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We  note  that  this  is  a  system  of  equations  for  two  holomorphic  functions 
on  a  strip  which  are  real  for  real  (the  bottom)  and  satisfy  a  pair  of 

coupled  nonlinear  time-dependent  boundary  conditions  along  the  top  of  the 
strip.  Kano  and  Nishida  [11]  used  essentially  the  system  (2.5),  along  with 
some  basic  facts  about  harmonically  conjugate  functions  on  a  strip,  to 
obtain  a  nonlinear  expression  for  the  t-derivatives  of  x  and  $  along 
n’  =  1  ,  for  which  a  solution  will  exist  to  the  Cauchy  problem  for  small 
times  (see  also  [18]). 

We  will  now  consider  a  particular  limiting  case  of  system  (2.7) 
corresponding  to  long  wavelength,  small  amplitude  waves  of  slow  time  variation 
and  will  obtain  the  Korteweg-deVries  equation  in  the  limit.  We  assume  that, 
as  |  C*  |  -*■  <x>,  -iO1  ■*"  0  and  z  "*■  1,  so  the  perturbations  from  the 
steady  flow  vanish  asymptotically.  The  limiting  case  is  given  by  the 
following  rescaling: 

Define  new  independent  variables 

(2.8)  C  =  aC  ;  n  =  n  ;  t  =  a^t  ,  where  a2  =  £, 
and  new  dependent  variables  x,y,G,A  by: 


ax(C,n,T,e) 

=  Re  ( z  ( ?  ’  ,  t )  )  -  £  ' 

2-  , 

a  y(f,,n,r ,f. ) 

• 

=  Im  (z  (C  ' ,  t)  -  o  • 

a2\  (r  ,r.,T,e) 

=  I'U'.t) 

( a  30  (f,n,T,<) 

=  ^  '  (’','  ,t) 
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Substituting  these  variables  into  the  system  (2.7)  gives  system  (2.3) 
above,  which  \;e  have  therefore  derived. 

B .  A  Formal  Solution  Procedure 

If  we  consider  the  system  (2.3)  and  assume  expansions  for  A,0,x,y 
of  the  form: 


A(r,n,T,e)  = 

*■* 

ej 

i  =  0 

e (C»h ,c,c)  = 

Of- 

9  ^  (C,n,T)  e  ^ 

j=o 

(2.10)  • 

00 

x(£,n,T,e)  = 

r« 

\ 

L 

x  ^  (Mh/T) 

j=o 

y(S,n»T,e)  = 

00 

I 

y  ^  (C»n»x) 

j=0 

then  (i) 

of  (2.3)  implies 

the 

system: 

(2.11)  • 

A<j)  +  0(j)  = 

£  n 

0  ; 

-x(j)  +  e^'15 

n  ^ 

=  0 

.  ;<3)  *  - 

0  ; 

-;(j)  + 
n 

=  0 

( 0 ) 

In  particular,  \  v  ~  0  ; 

n 

-  (0) 

x  =  0  .  From 

n 

the  1 

on  9 ,  y 

at  n  =  0  ,  this 

implies: 

»<°)  .  »«”(<;. T 

)  ? 

-(0)  '(0). 
x  =  x  ( 

(2.12)  1 

„<0)  ,1“',. 

-  -  F  l  ■  < 

r )  n 

"(0)  M0), 

;  y  =  x.  ( 

f ,  T  ) 
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SO  A  (1)  (£,n,T)  =  -  ~  (C,T)  n2  +  Q(1)  (C,T) 

x(1)  (E,n,T)  =  -  ~  (£,t)  n2  +  p(1)  (S,x) 

0(1)(f;,n,T)  =  n3-Q^(1)(C,T)  n 

-  (1)  ,r  ,  l  -  (0)  ..  .  3  (l)  „ 

y  (€ , n , r )  =  -  -  x  (C,t)  n  +  p.  (S,t)  n 

o  s 

Proceeding  inductively,  with  Q  ^  H  A^  ,  P^  H  x^  ,  we  have 


|  *  (-0^)  =  }.  — jSjjT-  (35)  2  «,t) 


(2.13) 


»(k)  ...  .  .  "  (-1)  jn2j  ,  3  .2j^(k-j) 

x  (f'n'T)  +  J  ■  (2j rr  ^  p  (C'T) 

3=0 


(i )  -'{-i) 

with  similar  expressions  for  0  ,  y  involving  odd  powers  of  n  and 

3/3?;. 

( i )  *  ( i ) 

Thus  if  A  ,  j  <  k,  are  known,  and  similarly  for  x  ,  there 

(k)  (k) 

are  two  unknown  functions  P  (f, ,  x )  ,  Q  (£,x)  which  arise  in  terms  of 
order  ck  and  higher.  Substituting  these  series  (2.13)  into  the  two 
boundary  conditions  at  n  =  1  ,  namely 


2<  '  -3 


k-  A  rlf  3/?  ?  a  a  3a  ^ 

y.  +  e  jcos(c  0)[cA^  +  c  (A tx ^  -  A^x^)  +  e  (0 t -  Qj-y^)] 


-"iP-  <-£\  *  -  »e?,  -  VE  ♦  V,”}  ■ 0 


along  n  =  1 


,  3/2,..  a  sin(c3/2fi)  -cXf  ~  ,  2, a  - 

o-  (  n)  -  (l+sxr)  - — -  +  e  <cy^  +  e  (x^yT 

c  ^ 


(Obtain,  setting  r  -  0, 


\V 


0  along  n  =  1 


( 2  .  1 4 ) 


from  rncli  <•  ;uat  ion. 


(°)  +  p(0)  =  0 


-10- 


Terms  of  order  r  in  the  boundary  conditions  at  n  =  1  are: 


(2.15) 


,<1)  +  2X(0),(0)  +  Ml)  .  3-(0)  rt(0) 
S  C  5  5  t 


=  0 


yc{1)  -  0(1) 


^(0)g  (0)  +  -(0)  =  Q 

?  T 


which  imply: 


1  (0)  (1)  (0)  (0)  1  _(0)  _(1)  ,_(0)  .  _(0) 

2  2^?  +  2c  +  22  2e  •  6  +  p?s  "  3Pcc  +  2t 


=  0 


P<1J  _  I  p<0>  _  I  Q(0)  +  Q(l)  +  p(0)  (0)  +  (0)  = 

so  that  P^'  +  is  known  in  terms  of  Q^,  P^  and  drops  out  upon 

subtracting  these  two  equations.  If  we  integrate  (2.14) ,  we  have 

+  p/°^  =  0  (by  our  boundary  conditions  as  |^|  -*•  °°)  ,  so  that  the  two 
boundary  conditions  for  order  e  imply: 

(2.16)  2Qt(0)  +  3Q(0)Q?(0)  +  3QC(0)  -  \  =  0 

which  is  a  form  of  the  KdV  equation. 

Remark The  formal  expansion  of  Friedrichs  and  Hyers  [9]  for  the  solitary 
wave  has  the  time-independent  form  of  (2.16)  as  the  equation  for  the  leading 
term. 

If  we  pick  any  satisfying  (2.16),  we  obtain  P^  by  integration, 

since  -  -  Q<(». 

The  order  e  terms  in  the  boundary  conditions  yield  two  equations 
of  the  form: 


(2.17) 


X<k)  +  2)(0,b(k-11  +  2X(k-1)Xr(0)  +  ;r(k)  -  Sy*-"  +X(k-1}  =  R 

f,  t  K  t,  6  t  k-: 


y<k)  -e(k)  -  ^o)0»-1)  -  xc(k"1)9(0)  + 


where  (and  later  R^  ,  S^)  depend  only  on  P^,  for  j  <  l . 
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(k)  (k) 

Thus  +  Q  =  i  and  again,  by  suotrustion, 

(2.18)  2A^k_1)  +  3(A(0)X(k_1))c  +  -  jA^J1*  =  Sk_2 

where  we  used  P^k  ^  +  Q^k  ^  =  R^_  ^  to  eliminate  x^k  ^  in  (2.18). 

Inductively,  we  find  a  formal  solution  using  the  power  series  (2.13)  by 

solving  (2.18)  for  A ^k  ^  and  then  obtaining  x ^k  ^  by  the  relation 
(k-1)  (k-1)  - 

pr3  +  Qa  ^-2* 

The  nontrivial  step  is  solving  (2.18),  the  linearized  KdV  equation 
with  inhomogeneous  terms.  For  water  wave  solutions  of  the  system  (2.3) 
which  behave  like  N-tuple  solitary  waves,  we  would  choose  for  A  ^  an 
N-soliton  solution  of  the  KdV  equation  and  then  solve  (2.18)  with  this 
X^,  seeking  solutions  with  the  appropriate  asymptotic  behavior. 

C.  Reformulation  as  a  Mapping 

We  may  formulate  the  equations  in  (2.3),  somewhat  artificially,  as 
components  of  a  mapping  F(e,0,y),  where  the  equations  (2.3)  correspond 
to  F(e,Q,y)  =0.  This  will  be  the  analogue  of  Beale's  approach  in  the 
stationary  case  [4].  We  begin  with  (2.3),  ignoring  the  fact  that  this 
is  the  rescaled  version  of  (2.8). 

Thus,  we  consider  functions  y(£,n,T)  and  0(£,n,i)  satisfying 
y  (5,0,t)  =  0,  0(S,O,t)  =  0  and  define: 


(2.19)  x  (f,  ,n  ,t)  = 


yTi(S,,n,r)  d£  *  j 


A(i,n,T)  =  - 


en(C’,n,x)  d£»  . 
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n 


The  vector  of  functions  F(e,0,y)  =  (F.^ ,  F^  ,F ^  jF^)  ■‘■s  as  foH°w'' 

(note:  F  ,  are  functions  of  £,,n,r  while  F^,  F^  depend  only  on 


F  =  c0  +0  ;  F  H  e y  +  y 

1  nn  2  fit,  nn 

_  .  2eX  a  ~3c  cXf  /  3/2A .  ,  ,  2,,  *  ,  A  . 

F,  =  X.e  +  y_e  +  e  Scos(e  0)[eX  +e  (X  x  -X  x  ) 

3£  C  (.  ttCCt 

(2.20)-  +  ‘WV/1 

*  sinlv2Ze)  '-eVe4(M5-VT-9TVeEV1} 

_  _  *■  .  3/2„ .  ,,  .  *  .  sin(c3//20)  ,  -eXf  A  ,  2  ,A  A 

F4  =  y^cos  (e  0)  -  (1+ex^)  - ^ -  +  e  jey^  +  e  (x^ 


xtV 


where  F^,  F^  are  evaluated  at  n  =  1. 

The  degeneracy  at  e  =  0  which  we  observed  in  the  formal  expansion 

( i )  a  ( i )  i+i 

above  (X  ,  x  were  only  determined  by  order  e  equations)  may  be 

removed  by  considering  a  modified  operator  F(c,0,y),  which  has  the  same 

roots  as  F  for  e  >  0. 

We  examine  the  kernel  of  F(O,0,y):  If  F(0,6,y)  =  0,  then  0^  =  0  ; 

y  =  0  ;  y  -  el  1  =  0  ;  y_  -  8  I  =  0.  Recalling  that  y|  =  9 1  =  0 

nn  i  n=i  C  n  n=i  n=o  1  ’-,=o 

this  gives  the  solutions: 


(2.21)  y  =  GU,t)*  n  ;  0  =  G  (£,i)  n 


where  G(f,,i)  is  arbitrary  except  for  the  asymptotic  condition: 

G(£,t)  ■*  0  as  |c|  -*■  00 .  To  study  the  range  of  F(O,0,y)  ,  suppose 

f1  f1 

F(O,0,y)  =  (a,6,r,s)  .  Then  since  0  —  0 |  ,  =  n*9  dn  =  n*a  dn, 

n  n=i  JQ  nn  JQ 

we  must  have  the  compatibility  condition: 


(2.22) 


f1 

|  n*a  dn  . 

1  o 


If  we  suppose  that  a,P,r  are  arbitrary  functions  which  vanish  as 
|f|  -»  °°  for  r  fixed,  and  define  s(r.,r)  by  (2.22),  i.e. 
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1 


s  =  r  + 


a  *o 


dn  then  F(O,0,y)  =  (a,0,r,s)  if 


(2.23) 


0  (£,n,f) 


y(C»n,t) 


o 

n 


(n-n')  a(£,n’,T)  dn*  +  A(S,T)n 

(n-nf)  B(S,n',T)  dn'  +  B(5,t)n 

l 


where  B^.(C/T)  -  A(C,t)  -  s(£,t) 


(1-n ' )  (B^_a)  dn' 


We  will  define  a  modified  operator  F  by  means  of  a  projection  Q  on 

the  range  space  which  incorporates  the  solvability  condition  (2.22).  Define 

1 

(2.24)  Q(o,B,r,s)  =  (0,0,0,  s  -  r  -  j  <x*n  dn) 

0 

Clearly  Q  is  a  projection  and  QF(O,0,y)  =  0.  Defining  F(e,0,y)  as: 


(2.25)  F  (e ,0 ,y) 


r  9?  <e » e_tY-l  +  (x-Q)  F  (e  ,0  ,y)  ,  t  >  0 
e 

QFe(0,6,y)  +  (I-Q)  F (0,0 ,y)  ,  e  =  0 


F  is  a  smooth  operator  with  the  same  roots  as  F  for  e  >  0,  but  the 

additional  condition  QF£  =  0  leads  to  more  regular  behavior  at  e  =  0. 

In  fact,  the  added  condition  is  precisely  the  KdV  equation  i.e.  suppose 

F  (0,0 ,y)  -  0.  Then  QFc(O,0,y)  =  0  means  precisely: 

1 

(2.26)  yT  -  ox.  -  XT  +  3yc  -  2\\(.  -  j  6^-n'  dn'  =0  for  n  =  1. 

0 

Since  QF  -  0  at  e  =  0,  (I-Q)F (0,9 ,y)  =  0  implies  0  =  G?(C,t), 
y  -  n(r,T)ii.  Using  the  relation  A  =  -  I  0^  d?  =  -G(£,t)  in  (2.26) 

•  — oo 

above,  we  obtain  an  equation  for  A  ,  namely: 


(2.27) 


2A 


-  3AA  -  3A 

t. 


+  3  ° 
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where  X  \(£,x)  ,  y  -  -X(£,x)n  ,  0  =  -X^(^,i)n.  This  is  the  KdV  equation 

(2.16)  above. 

In  the  case  of  the  solitary  wave,  Beale  was  able  to  find  an  appropriate 
scale  of  Banach  spaces  for  the  domain  and  range  of  F.  He  then  proved  that 
in  a  neighborhood  of  e  =  0  ,  9  =  (corresponding  to  the  1-soliton  KdV 

solution),  the  Fredhet  derivative  of  F  was  invertible  in  the  generalized 
sense  of  Nash-Moser.  This  implied  the  existence  of  a  nontrivial  solution 
for  e  >  0  sufficiently  small  by  means  of  a  generalized  implicit  function 
theorem  (Zehnder  [25]). 

For  the  problem  considered  here,  we  have  not  as  yet  specified  the 
proper  function  spaces.  Continuing  on  a  formal  level,  we  consider  the 
linearized  operator  at  j  =  0,  9  =  0Q  ,  y  =  yQ  ,  where  F(C,0Q,yo)  =  0. 

A 

We  shall  define  another  projection,  P,  mapping  (0,y)  onto  the 
kernel  of  F(0,* ,*)  by: 


(2.28)  p  (9  ,y)  =  (y^(C,l,T)-n,  y(C,l,x)‘n) 


We  note:  P(9,y)  =  (0,0)  if  and  only  if  y(£,l,x)  =  o.  Using  the  pro¬ 
jections  P  and  Q  ,  we  may  regard  the  linearized  operator  dF  Q  ~  . 

(0,90,y0) 

as  a  matrix  of  operators: 


(2.29) 


dF  = 


Q  dF|p  (I-Q)  dF|p 

QdF|x_r  (I-Q)dF|l_p 


Since  (I-Q)  dF |  1  0  ,  invertibility  of  dF  reduces  to  invertibility  of 

the  'diagonal  elements'  Q  dF]p  and  (I-Q)  dF|x  p.  We  remark  that 

(I-Q)  F I  is  a  linear  operator,  whose  inverse  we  may  compute  explicitly. 

A  ”  i  _ 


The  equation  (I-2)F|I_p  (0,0, Y) 
solution : 


(a,B,r,r+ 


na  in) 
o 


has  the  unique 
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(2.30)  0(£,n,x)  i 

Y(£,n 


(n-n*)  a(C,n'»T)  dn'  -  n[r(E,r)  + 


a(6,o',T)  dn’) 


,t)  =  [  (n-n' 

J  o 


)  ea,n',T)  -  m 


(l-n')B(5,n\T)dn'] 


o 


which  is  obtained  by  adding  the  condition  Y(£,1,t)  =  0  to  the  solutions 
given  previously  in  (2.2  3). 

The  essential  difficulty  with  this  inversion  is  that  the  Banach  spaces 
for  0fY  may  include  ?  and  t  derivatives  in  their  norms;  for  this 
reason,  Beale  uses  a  generalized  implicit  function  theorem  in  the  station¬ 
ary  case  (shrinking  the  domains  of  analyticity  in  order  to  control  deriv¬ 
atives)  . 

The  second  diagonal  element,  Q  dF|p  ,  is  invertible  if  we  can  solve: 
(2.31)  -  2Hfc  -  3H?  -  3  ( X  »H)  +  j  H  =  J(£,x) 

the  inhomogeneous  linearized  form  of  the  KdV  equation  (see  equation  (2.18) 
above) . 

For  the  remainder  of  this  paper,  we  consider  the  linearized  KdV 
equation.  As  we  have  seen,  it  arises  in  the  study  of  small  amplitude, 
long  wavelength,  slow  time  variations  of  a  steady  flow  of  a  perfect  fluid 
over  a  flat  bottom  with  Froude  number  near  1.  If  we  seek  solutions  describ¬ 
ing  a  'nonlinear  superposition'  of  N  solitary  waves  of  distinct  speeds, 
the  first  approximant  will  be  an  N-soliton  solution  of  the  KdV  equation 
ind  the  higher  order  corrections  will  satisfy  the  inhomogeneous  form  of  the 
linearized  KdV  equation  (linearized  about  the  N-soliton) . 

We  shall  consider  the  Cauchy  problem  for  the  linearized  KdV  equation. 
By  Duhamel's  principle,  this  amounts  to  solving  the  inhomogeneous  equation. 
By  the  change  of  variables, 
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(2.32) 


X  -  9T  =  £ 


-  6T  =  T 

q(X,T)  =  |  X(£,r) 
we  obtain  the  usual  form  of  the  KdV  equation 

(2'33)  qT  +  qXXX  -  6q<?X  =  ° 

We  note  that  the  T-independent  solution  of  the  KdV  equation  is  a  function 
of  X  -  9T  ;  this  gives  the  one  soliton  with  speed  9,  which  explicitly  is 
X(0)  =  -3sech2(-j  £)»  the  first  order  term  on  the  expansion  of  Friedrichs 
and  Hyers. 

In  the  remaining  sections,  we  shall  use  the  letters  x,y,t,u,v  etc. 
for  meanings  other  than  those  of  the  above  section.  Since  these  different 
meanings  occur  in  separate  places,  this  should  cause  no  confusion  for  the 
reader. 
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3.  SOME  RESULTS  ON  THE  CAUCHY  PROBLEM  FOR  THE  LINEARIZED  KdV  EQUATION 


In  this  section,  we  summarize  the  results  of  [19]  regarding  the 
Cauchy  problem: 


(3.1) 


where  q(x,t) 


u 

XXX 


-  6<qu)x 


0 


[  u(x,0)  =  4>  (x) 
satisfies  the  KdV  equation 


(3.2) 


+  q 

XXX 


-  6qq  =  0 
x 


By  Duhamel's  principle,  the  inhomogeneous  form  of  (3.1)  is  solvable  if  the 
Cauchy  problem  is. 

In  [19],  an  explicit  formula  for  the  solution  of  problem  (3.1)  is 
given,  using  certain  functions  arising  from  the  SchrSdinger  equation 

(3.3)  -fM(x,k,t)  +  q(x,t)  f  (x,k,t)  =k2f(x,k,t) 


where  the  potential  q(x,t)  satisfies: 


(3.4) 


2 

(l+x  )|q(x,t)|dx  <  00  for  every 


fixed. 


The  fundamental  discovery  of  Gardner,  Greene,  Kruskal,  and  Miura  [10], 
later  formulated  abstractly  by  Lax  [15],  is  that  if  q(x,t)  evolves 
according  to  the  KdV  equation  (3.2),  the  spectrum  of  the  Schrbdinger  equation 
(3.3)  is  fixed  and  the  associated  scattering  data  evolves  in  a  simple  way. 

We  shall  use  this  information  below,  but  first  introduce  some  notation 
and  basic  facts  about  the  scattering  theory  for  (3.3).  This  information 
(and  much  more)  may  be  found  in  [7] . 

Let  f^(x,k,t)  denote  the  Jost  solutions  of  (3.3) 

1 

,  ,  ,  „  lkx  +  4ik  t  ,  ^  ...  . 

i.e.  t+(x,k,t)  a.  e  as  x  +°°  ,  t  fixed 

f  (x,k,t)  'v  e  :''kX  t  as  x  ->  -00  ,  t  fixed 
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and  both  satisfy  (3.3).  We  define  the  transmission  coefficient,  T(k,t), 
in  terms  of  the  Wronskan  of  f+,  f  as  follows: 

(3-5)  TTkTtT  =  2lk  [f+(x'k't)'  f_ (x,k,t) ] 

f|(x,k,t)  f_(x,k,t)  -  fMx,k,t)  f+(x,k,t) 

=  __ 


3  3 

(We  shall  always  use  the  notation:  *  =  —  ,  *  =  —  .)  It  is  not  hard  to 

o  X  dk 

show  that  T (k , t)  =  T(k)  is  independent  of  t  and  that  under  the  normal¬ 
ization  of  f+,  f_,  T(k)  is  meromorphic  in  the  upper  half-plane  Im  k  >  0 

2 

with  poles  at  k  =  if? ^ ,  j  =  1, . . .  ,N  where  each  energy  is  a  bound 

state  energy  for  (3.3).  N  is  finite  by  a  classical  estimate  involving 
(1+| x| ) | q(x) 1 dx  <  °°.  T(k)  is  also  continuous  and  non-zero  for  real 

k  f  0  .  For  notational  ease,  we  introduce  for  j  =  1,...,N  the  following 
pair  of  functions: 


(3.6)  F  j  (x,  t)  =  f+  (x,i8  j  ,t)  ;  G^  (x,t)  =  c^  f+  (x, iB ^  , t)  •  g..(x,t) 


where  g.(x,t)  =  [f  (x,k,t)  - 

3  l  dk  - 


f_(x,iB . ,t) 

-z~. - .  rP-  f.(x,k,t)]i.  and  c. 

f+(x,iB  ,t)  +  ] K=iB ^  3 


F((x,0)  Gj(x,0)  dx  =  1  for  j  =  1,...,N. 


3  1  3 

r 

is  chosen  so  that 
The  principal  result  of  [’.9]  is  the  following: 

Theorem  3.1  Suppose  q(x,t)  satisfies  (3.4).  If  <p  (x)  is  continuous 
and  integrable,  the  solution  of  (3.1)  (in  the  sense  of  distributions)  is 
given  by: 


(3.7)  u(x,t)  = 


_d]L_  T2(k)  j 
4:rik  1  '  1 


co 

j  |  [f^(x,k,t)ff (y,k. 


0) 


-  f2_  (x,k,t)  f^(y,k,0)  ]  <J>(y)  dy 


N 

l 

3*1 


[f!  (x,t)Gj  (y,0)  -  c/  (x,t)F.  (y,0)]  0>(y)  dy 
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\ 


For  a  proof  of  this  theorem,  see  [19] . 

Remark :  It  is  known  ([10],  Theorem  3.6  or  [19])  that  the  functions 
2 

(f+) '  (x,k,t) ,  F!(x,t),  G! (x,t)  all  satisfy  the  linearized  KdV  equation 
(3.1).  The  formula  for  u(x,t)  resembles  the  Fourier  decomposition  of 
< p  (x) ,  where  the  derivatives  of  the  squared  eigenfunctions  replace  the  usual 
exponentials  and  the  presence  of  a  non-zero  potential  q(x,t)  can  lead  to 
the  discrete  terms  F*(x,t),  G!(x,t).  In  fact,  when  q(x,t)  H  o,  (3.7) 
reduces  to  the  usual  Fourier  transform  solution  of  the  Cauchy  problem: 


namely 

CD  OO 

3  (  \ 

2ikx+8ik  t  J  -2iky  ,,  .  , 

S  e  ()>(x)  dy  ■ 
l  > 

—00 

3 

,  _  ,i  .  _  -  ±(ikx+4ik  t) 

since  for  q  =  0  ,  T (k)  =  1  and  f+(x,k,t)  =  e 

Noting  that  the  solution  u(x,t)  given  by  (3.7)  consists  of  two 
pieces —  a  discrete  sum  and  an  integral,  we  analyze  them  separately.  The 
sum  corresponds  to  variations  in  the  soliton  part  of  the  function  q(x,t) 
and  decomposes  into  travelling  waves  with  positive  velocities  as  t 
becomes  large.  For  the  water  wave  problem  of  Section  2,  these  terms  are 
of  considerable  interest.  The  k-space  integral  part  of  (3.7)  forms  a 
dispersive  wave  train  and  will  be  seen  to  behave  like  the  solution  v(x,t) 
of  the  Airy  equation  (3.8).  In  particular,  for  initial  data  which  is 
somewhat  smoother  and  more  rapidly  decaying  than  was  assumed  in  Theorem  3.1 
above,  we  show  that  this  part  of  the  solution  u(x,t)  is  smoother  for 
t  >  0  but,  as  x  -*■  -  °°,  it  decays  less  rapidly.  We  present  these  results 
for  linearizations  about  N-soliton  solutions  of  the  KdV  equation.  Similar 
analysis  applies  for  a  more  general  class  of  KdV  solutions  satisfying 


v(x,t)  =  — 

TT 
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* 


(3.4) ;  we  omit  such  a  discussion  for  the  sake  of  brevity  and  restrict  our 
attention  to  the  N-soliton  case. 

Slower  decay  as  x  ■>  -»  for  t  >  0  occurs  because  of  the  dispersive 

d  2 

nature  of  the  oscillating  solutions  —  (f^(x,k,t))  of  the  linearized  KdV 

equation  (see  [24]  for  a  general  discussion  of  dispersive  waves) .  In  part- 

2 

lcular,  the  asymptotic  behavior  of  f+(x,k,t)  ,  as  x  -*•  ±®>  respectively, 
is  given  by  the  exponentials  (x,k,t) ^  w^ere  we  aefine: 

(3.10)  0(x,k,t)  =  2kx  +  8k ^t. 

2 

These  waves  propogate  with  a  negacive  velocity  -4k  so  that  waves  with 
large  wave  numbers  contribute  to  the  solution  near  x  =  -ro  almost 
instantaneously. 

i  i0 

The  same  expontentials ,  e  ,  form  the  solution  of  the  linearized 
equation  for  q  =  0  ,  (equation  (3.9)  above),  namely 


V  +  V 
t  XXX 


0  , 


as  is  seen  by  Fourier  transform,  and  arise  in  the  asymptotics  of 
2 

f+(x,k,t),  which,  by  the  trace  formula  of  Deif t-Trubowi tz  [7],  lead  to  a 
solution  of  the  full  KdV  equation 


(3.11) 


q.  +  q  -  6qq  =  0. 

t  XXX  X 


(In  [7],  q  (x , 0)  is  written  as  an  integral  over  the  real  line  in  k: 

OO 

•  .  n 

(3.12)  q  (x , 0)  =  — ■  .  (k,0)f^(x,k,0)dk  +  T  a . f2 (x , i6  .  , 0) . 

'■  +  j  =  ]  '  ’  ' 

— .oo  J 

An  approach  to  tie  KdV  equation  itself  using  (3.12)  will  appear  in  a 
subsequent  paper  by  the  author.)  The  smoothness  and  decay  properties  of 
the  solution  of  the  Cauchy  problem  for  the  KdV  equation  were  analyzed  by 
Tanaka  [23]  and  later  Cohen  Murray  [5]  using  Faddeev-Marchenko  i nverse 
■‘•uttering  theory  rather  than  the  then-unknown  trace  formula  (3.12); 
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asymptotic  analysis  of  the  KdV  equation  also  appeared  in  [1,2],  where  the 
more  delicate  regions  —  =  0(1)  as  t  -*•  +  °°  were  also  discussed  in  the 
absence  of  solitons. 

Our  analysis  for  the  linearized  KdV  equation  proceeds  in  direct 

analogy  with  the  equation  (3.8);  the  chief  difference  is  the  presence  of 
2 

the  factors  m+ (x,k,t)  multiplying  the  exponentials  and  their  derivatives, 
which  must  be  considered  in  all  arguments.  The  techniques  used  will  be 
primarily  integration  by  parts  and  stationary  phase  analysis.  In  the 
limits  we  consider,  the  stationary  phase  points  tend  to  ±°°  ,  which 
complicates  matters  slightly.  As  in  [5] ,  we  will  work  in  a  shrinking 
neighborhood  of  the  stationary  phase  points,  whose  size  is  proportional  to 
a  small  negative  power  of  ]x|.  This  variation  of  the  usual  stationary 
phase  argument  [8]  is  used  to  control  the  error  terms  arising  at  the 
stationary  phase  points.  The  smoothness  argument  relies  on  the  observation 
[3]  that  for  t  >  0,  we  may  rewrite  the  x-derivative  of  0  in  terms  of 
the  k-deri vati ves  of  0  as  follows: 


(  3.13) 


(n  )2  -  4k2  =  . 

x  6t  3t 


We  will  use  this  to  re-express  u  as  a  function  which  is  smoother  than 

;  mi  otherwi  so  appear  to  be. 

X  '< 

re  >ults;  are  summarized  in  the  following  theorem: 

THi-.i  ■KKM  > .  2  Assume  1  (x)  ,  the  initial  data  for  the  linearized  KdV 

!•.;  i  ••  :  has  four  continuous  derivatives  and  that,  for  some  fixed 


(1  !  x 


Mx)  G  L  for  0  <  a  £  4. 


'  when  q(x,t)  is  an  N-soliton  solution  of  the  KdV  equation, 
mod  in  (3.7)  above,  has  the  following  properties: 
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(3.15) 


i 


(i)  u(x,t)  is  a  classical  solution  of  (3.1)  for  t  >  0 
with  u  (x , 0)  =  <f>  (x) 

3T  S 

(ii)  3  3  u(x,t)  is  continuous  for  t  >  0  for  all  non-negative 

integers  r,s  satisfying  3r  +  s  <  21  +  2 

(iii)  lim^J  u(x,t)  •x*' [  =0  t  >  0  fixed 

(iv)  | x |  |u(x,t)|  is  bounded  as  a  function  of  x  for  t  >  0 

(even  as  x  -*•  -«) 

(v)  |u(ct+6,t)  |t  is  bounded  for  c  <  0  as  t  -*■  +“>,  5  fixed 


The  proof  of  Theorem  3.2  is  given  in  the  three  lemmas  below,  in  which 
the  smoothness  and  the  limiting  behavior  are  discussed  separately.  First, 
we  present  some  facts  concerning  the  Jost  functions  f+(x,k,t)  in  the 
N-soliton  case,  where  we  choose  the  phases  of  the  waves  so  that  q(-x,-t)  = 
q  (x , t) .  (Recalling  the  scaling  done  in  Section  2,  this  normalization  is 
reasonable. ) 

The  explicit  form  of  the  N-soliton  leads  to  an  algebraic  expression 
for  the  Jost  functions  (see  also  [6]).  In  the  proof  of  Theorem  3.2,  we 
shall  exploit  certain  properties  of  these  functions,  which  we  state  here 
and  prove  in  the  Appendix.  Define,  for  j  =  1,...,N, 


(3.16) 


cosh (8 .  £ . ) , 
3  3 

r  .  x  -  48 2 1  ;  i’.  =  - 
3  3  3 

sinh  (P_.  £  ^  , 


j  odd 

j  even 


and  consider  the  N  ;<  N  Wronskian  determinant  (in  x  ): 


fixed 
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(3.17)  w(x,t)  =  det 


1 

(1) 


<P 


4- 


2 

(1) 


ip  (N-l)  (};  (N-l) 

*>  1  2 


W*"'V 


4> 


'P 


N 

(1) 

N 


\P  (N-l) 
N 


In  the  Appendix,  we  show  w(x,t)  >  0.  The  N-soliton  solution  of  the  KdV 
equation  is  given  by: 

d2 

(3.18)  q  (x , t)  =  -2  — -  log  w(x,t) 

dx  , 


while  the  eigenfunctions  f+(x,k,t)  are  given  by  the  ratios: 


w. 


(3.19)  f+(x,k,t)  = 


N+l 


^  1 '  ^ 2 '  "  ‘  '  '^N,e 


±i  (kx+4k  t) 


w(x, t)  ||  ± (ik  -  0 . ) 

j=l  3 


where  is  the  (N+l)* (N+l)  Wronskian  determinant.  Writing 

ff(x,k,t)  =  m+ (x ,k , t) ^  ,  we  deduce  the  following  properties  of 
the  factors  m+(x,k,t)  from  (3.19): 


(3.20)  (i)  m+(x,k,t)  are  rational  functions  of  k.  Their 

denominations  and  numerators  are  polynomials  of  degree 

N  in  k  ;  both  denominators  are  in  fact  precisely 
N 

"|  f (k  +  iS  ■ )  while  the  numerators  are  polynomials 

j=l  3 

N  +  ,  N-l  +  + 

k  +  A  (x,t)k  +  ...  +  A  (x,t)k  +  A  (x,t)  where 
1  N-l  N 

4- 

each  coefficient  A  (x, t)  is  a  rational  function  of 

f  B  1 

[e  3  3j  which  is  bounded.  The  denominator  of  each 

4- 

A  _(x,t)  is  w(x,t),  which  we  show  in  the  appendix  is 
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L  e.B.C. 

a  sum  of  terms  eJ  *  over  all  possible  choices 

e =  ±1  with  positive  coefficients  for  each  term. 


(ii)  —  m+(x,k,t)  is  a  rational  function  of  k  which  decays 


like  k 


as  |  k  |  -*■  °°. 
N 


(iii) 


_d_ 

dx 


m+(x,k,t)  = 


N  (  T 

£  fe]  (A,±(x,t)]-kN-J 


£=1 


TT<*  +  i6.) 

j=l  1 


so  all  x-derivatives  of  m+  decay  like  |k 
!k|  ■>  “. 


-1 


N  k+iB . 

(iv)  We  also  have:  T(k)  =  I  [  : — 

.  '  k-iB . 

3  =  1  3 

(3.21)  F j  *  (x , t)  and  G^'(x,t)  are  real  analytic  in  x,t  and  for  fixed  t 
they  decay  exponentially  fast  as  |  x  |  -+  »  (see  Appendix). 

By  (3.21),  all  the  smoothness  and  decay  properties  of  Theorem  3.2 
are  satisfied  by  F^'(x,t)  and  G^'(x,t) .  Therefore,  we  consider  the 
function  H(x,t)  given  by: 


(3.22) 


Q  (x,  t) 


f  ^  (k)  d  f  ?  2  2  2 

I  4irik  dx  jf+  (x,k,t)  f_  (y,k,0)  -  f_  (x,k,t)  f^  (y,k,0) 


—00  — oo 


(y) dy  dk. 


Note  that  the  integrand  is  continuous,  even  at  k  =  0  (since  f  (x,0,t), 
f_  (x,0,t)  are  linearly  dependent).  Formula  (3.22)  suggests  the  following 
definitions: 


(3.23) 


$+(k)  H 


f+  (y,k,0)<f>  (y)  dy 


00 

m2  (y,k,0)e±2lky<j,  (y)  dy 

—00 
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Wo  will  analyze  ij>t(k)  just  as  in  the  usual  Fourier  transform  case, 

using  (3,20)  to  control  the  extra  terms.  Thus  we  shall  see  that  u(x,t), 

given  by  (3.22',  and  v(x,t),  the  solution  to  the  linearized  problem  for 

q  E  0  given  by  (3.9),  behave  quite  similarly. 

The  first  part  of  Theorem  3.2  is  contained  in  the  following  lemma. 

.  r  j  ■>  a  . 

Lemma  3 . 3  If  (1  +  j j  ^ )  I— j  <p  (x)  6  L  for  all  0  <  a  <  4  where  l  >  4 

2*  S  ^ 

is  fixed,  then  the  functions  3^  u(x,t),  where  u(x,t)  is  given  by 

(3.17)  are  continuous  for  all  non-negative  integers  r,s  satisfying 

3  r  +  s  <  2  i  +  2 . 

Proof  of  Lemma  3.3.  The  idea  of  the  proof  is  a  follows:  we  show  that 

5((k)  decay  rapidly  enough  as  |kj->  ®  that  we  can  dif ferentiate  (4.9) 

twice  with  respect  to  x  and  still  have  a  convergent  integral.  Then, 

2 

using  (3.13)  to  eliminate  the  -4k  factor  arising  from  the  exponentials 
and  the  estimate  (3.20)  (iii)  to  control  derivatives  of  m  (x,k,t),  we  show 
that  the  integral  for  can  be  differentiated  twice.  Repeating  this 


argument,  we  obtain  the  desired  result. 


Step  1.  We  show  that 
for  0  '■  '  s  i. 


dk 


-4. 


<jj^k)  exists  and  is  0  ( | k ]  )  as 


Proof  of  Step  1.  If  we  integrate  (3.23)  by  parts,  we  have: 

1 


(3.24) 


a  ,  (k)  - 


l  2ik 


d  ,  2 ,  ,  ...  ±2iky , 

(m+  (y,k,0)<*>  (y)  )e  dy. 


dy 


The  integral  is  absolutely  convergent  by  our  assumptions  on  $  and  the 
properties  of  m,  (y,k,0)  listed  above.  In  fact,  we  may  integrate  by 
part::  four  times  with  respect  to  y,  obtaining 


(3.2I>)  (k) 


d  I  4  .  2 


(nT  (y,k,0)<j>  (y)  )c±2lkydy 


( 2  ik ) 

—to 

and  the  integral  is  still  absolutely  convergent.  Thus  $  (k)  is 

t M j  k |  S  as  j  kj  *■  on  . 

.dnee  (1  +  |  x  | ‘")  if)  (x )  £  l\  the  k -derivatives  of  $  (k)  of  order  less 
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than  or  equal  to  l  all  exist.  Integrating  these  expressions  by  parts 
four  times,  we  find 


(3.26) 


Y 

$±(k)  =  0  ( |  k  | 


-4. 


as  |  k  | 


00  for  0  S  y  -  ^ » 


which  completes  Step  1. 

Step  2.  (The  smoothness  for  t  >  0) .  Writing 

00 


(3.27)  u(x,t)  H 


dk 

4ffik 


T2(k)^  [m2(x,k,t)el9$_(k)) 


—03 


(m2 (x,k,t)e  10$+ (k) ) ] 


we  have  a  continuous  integrand  which  decays  like  |k|  ^  as  |k|  -*■  “>. 
Therefore  we  may  differentiate  twice  with  respect  to  x  and  still  have 
a  convergent  integral.  This  yields: 


(3.28) 


UxX(x,t) 


dk 

4irik 


T2(k) 


It)  fv(x'k' 


t)el9$_(k)  -  m2(x,k,t)e  l9$+(k) 


Since  0^  =  2k  and  m+* (x,k,t)  decays  like  |k|  ^ ,  the  terms  on  the 
integrand  in  which  m+f  (x,k,t)  or  a  higher  derivative  appears  all  have 
decay  like  |k|  4  or  faster;  these  terms  can  therefore  be  differentiated 
twice  more  with  respect  to  x.  The  remaining  terms,  in  which  the  exponential 
is  differentiated  three  times,  are: 


(3.29) 


^  (~4k2) 

2tt 


jm2  (x,k. 


t)e'''9^  (k)  +  m2(x,k,t)e  (k)  dk 


T2(k) 

2tt 


~  il]  f^x,k,t)el9*~_(k) 


+  m2(x,k,t)e  i9$+(k) 


-7-  im2  (x,k,t)el9<f»  (k)  +m2(x,k,t)e  l9J  (k) 
2tt  3t  (  +  -  -  + 


dk 

by  (3.13) 
i  dk 


00  00 


12Tiit 


d  ,_2  , 


- 

[dk 


|~(T‘(k)m:(x,k,t)*  (k)) 
[_dx  +  — 


(T2  (k)m2 (x,k, t) (k) ) 


'‘I 


dk 
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where  we  integrated  by  parts  in  the  second  term.  Each  of  these  terms 
has  continuous  k-integrands  which  decay  like  |k|  ^  or  better;  hence 
they  are  also  twice  differentiable  with  respect  to  x.  Therefore  u(x,t) 
has  in  fact  four  continuous  x-deri vatives  for  t  >  0.  Repeating  this 
argument  iteratively,  we  obtain  u(x,t)  has  21  +  2  continuous  x-derivatives 


(since  we  can  only  bound 
argument  l  times) . 


dk 


$+(k)  for  Q  <,  y  <  i ,  we  may  repeat  the 


To  handle  t-derivati ves ,  we  note  that  differentiating  directly  in  t 
brings  down  a  factor  9  ^  H  8k3,  which  does  not  a  priori  lead  to  a 
convergent  integral.  However,  multiplication  by  8k3  may  be  expressed 

,  .  ,3 

in  the  sense  of  distributions  as 

'a'3 


3x 


plus  convergent  integrals.  Since 


3x 


u(x,t)  is  continuous,  so  is  u^_.  The  equation  (*)  then  gives 


higher  regularity  and  the  desired  result.  This  proves  Lemma  3.3. 
The  decay  as  x  -*•  +“,  t  >  0  fixed  and  finite  is  given  by: 


Lemma  3.4.  For  t  >  0  finite,  fixed,  ]u(x,t)  x^  ]  •+■  0  as  x  -*■  +«■>. 

Proof.  Once  again,  we  need  only  consider  u(x,t)  since  F^’(x,t),  G_.'(x,t) 

2 

decay  exponentially.  For  x  >  0,  t  >  0  we  note  that  0^  =  ~x  +  24k  t  >  0 
and  in  fact: 


(3.30) 


°kk 

48kt 

< 

f 

X 

ek 

2x  +  24kV 

12tJ 

-1/2 


'dote  also  0,  ,  ,  =  48t,  which  is  bounded, 

kkk 


(3.31) 


u(x,t) 


00 

p(x,k,t)ei9dk 

. 

—  OO 


by  defining 


Write 
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(3.32)  P(x.k,t)e10  =  ~ ~  |n2  (x ,k , t) el6J _  (k ) 


+  T2 (-k)  d  r  2 ,  ,  _  i6r  ,  ,  .  , 

4>,(-k) 

47rik  dx  +  _ 


-4 


We  note  p(x,k,t)  is  continuous  in  k  and  decays  like  j  k  [  as  |k|  -*■ 


as  does  l^j  p(x,k,t)  for  0  <  y  <  i. 

Integrating  (3.31)  by  parts  2.  times  in  k,  we  have 


(3.33)  u(x,t)  =  (i) 


JL  J_ 

3k  6, 


P  (x,k,t) 


e19  dk. 


Using  (3.30)  and  the  obvious  bound 
estimate,  for  t  >  0  fixed, 


<  —  for  x  >  0  we  obtain  an 


(3.34)  |u(x,t)|  <  (2x)  £C(t)  for  x  SM>  0 

1/2 

where  C(t)  is  polynomial  in  t  of  degree  at  most  2.-1. 

Moreover,  since  the  integrand  in  (3.33)  is  integrable,  by  a  simple 
modification  of  the  usual  Riemann-Lebesque  lemma  (namely,  pick  <  with 


then  approximate  by  a  smooth  function  and  integrate  by  parts) , 


_oo  « 

SL~ 

we  can  show  x  u(x,t)  -*■  0  as  x  •+  +“>  for  t  >  0  fixed,  which  proves 
lemma  3.4. 

Finally  we  discuss  decay  as  x  -*■  -°°  for  t  >  0  fixed: 

i  1 9/4 

Lemma  3.5.  As  x  -*  for  t  >  0  fixed,  u(x,t)*|x|  remains  bounded. 


Proof  of  Lemma  3.5. 

0,  =  2x  +  24k2t  =  0 
k 

Let  k  =  cik.  Then 

(3.35)  u(x,t)  = 


Define 
for  k 


-x 


12t 


1/2 


For 


±a.  As  x  -<*>, 


x  <  0,  t  > 
(for  fixed 


OO 

p  (XjkjtJe^dk 

— OO 


a 


OO 

,  iX0 

p  (x,aK  ,  t)  e  d« 

— OO 


0, 

t  >  0)  , 


a  -*■  +°°. 


3/2 


where  \  S  2- 


(12t) 


1/2 


and  9 (< ) 


-<  +  <  /3,  so  that  0(k)  =0  for  k  =  +1. 
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comes  from  the  terms 


1+e 


o  (x,a,t) 


iX[-K+<  / 3 ] 
e  dK  and 


-1+e 

r 


,  . ,  iX0 

p  (x,a,t)e  dK 


where  the  k  value  is  frozen  at  k  =  ±1  in  the  function  p.  The  extra 

term  [p(x,«K,t)  -  p  (x,a , t)  ] e^^dK  is  of  lower  order  for  large  X  by 

bounds  on  the  derivative.  In  the  case  considered  here,  (*->•“>  so  this 

error  term  may  become  large.  To  counteract  this,  as  in  [5] ,  we  consider 

a  very  small  interval  about  the  stationary  phase  points  k  =  ±1,  of  order 
i  -u 


x 


for  instance.  We  estimate  as  follows: 


(3. 36) 


1+r. 


p  (x,uk, t)e  dK  =  a 


l+£ 


p  (x,a,t)e'i'',l6dK 


1+e 


(p  (x,a< , t)  -  p (x,a,t) )elX0dK 


1+e 


=  up (x,a,t) 


e  dK 


1+e 


J  2d  ,  -  , 

a  *dk  P  {x'aic't) 


,  .  ixe 

(K-l)e  dk 


for  some  k  £  [1,  1+e] 


The  first  term  is  estimated  as  in  the  usual  method  of  stationary  phase; 

the  second  term  leads  after  an  integration  by  parts  to  an  estimate  of  the 

form  C-  *">  \  since  ^-P  decays  like  ]k|  With  e  =  ]x|  U,  0  <  \j  < 

this  term  is  of  order  |x|  U  as  jx|  ->  oo.  The  first  term  decays  like 

•  >  :  1-9/4  ,  1 

u  •  ■  ‘Hjx|  )  as  [ x  |  oo  and  is  the  leading  term. 

Similar  estimates  hold  on  the  intervals  [l-c,l],  f — 1— e , — 1 ] ,  [-1,-1+c] 


1/2, 


On  the  interval  [l+e,°°),  we  estimate  as  follows: 


(3.37) 


,  iX0,  -a  e 

p(x,aie,t)e  die  =  - 


iX[-2/3+e2+e3/3] 


iXe  (2+e) 


p  (x,a (1+e) ,t) 


1+G 


a 

iX 


00 

— 

f  ' 

~ 

d 

p  (x,aie,t) 

dk 

2  _  ! 

1+e 

_ 

K  1  J 

iX§ 

e  die. 


which  leads  to  a  bound  of  the  form: 


(3.38) 


ixe . 

p  (x,aie,t)e  die 


1+e 


5  (X'VVV  C, 


,-l  -2  -1  „  ^  ,-l  -3  -2  „ 

+  X  a  e  *C2  +  X  a  e  *c 

=  0  ( [x [  3+2u)  since  0  <  u  <  1/2 


A  similar  estimate  holds  on  the  interval  (-°°,-l-e]. 
Finally,  on  the  interval  [-l+e,l-e],  we  have: 


1-e 


(3.39) 


p (x,aie,t)e 


iX0d|<  _  a«p  (x,aic,t)e 


ixe 


-1+e 


iX  (ie  -  1) 


ie=l-e 


ie=-l+e 


a 

iX 


1-e 

r 

“T 

d 

p (x,ax,t)> 

die 

2  -  1 

1+e 

_ 

K  J 

Integrating  by  parts  three  more  times,  we  obtain: 
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1-e 


(3.40)  a 


/  iX9. 

p  (x,ax , t) e  dx 


-1+e 


3 

l 

1=0 


-1 


[iX (<  -1) 


_d_ 

dx 


(x,ax , t) 
iX(K2  -  1) 


ixe 


x=l-e 


X=-l+G 


.-e 

4 

a 

d  1 

p  (x,ax,t) 

dx  ...  2  , . 
iX  (x  -1) 

-1+e  *■ 

.V  * 

iX0 

e  die . 


The  boundary  terms  lead  to  estimates,  for  &  =  0,1, 2, 3  of  the  form: 


,-l  -3  -1  . -2  -3  -3  , 

X  a  e  ;  X  a  e  (1  +  ac) ; 


,  -3  -3  -5„  ,  2  2,  .  -4  -3  -7  ,,  J  ^  2  2  ^  3  3. 

X  a  e  (1  +  ac  +  a  g  )  ;  X  a  e  (1+ae+ae  +  a  g  ) 


respectively. 

Since  a  is  0(|x|1//2),  X  is  0 ( | x |  3/^2 )  and  e  =  |x|  U,  the 

inequality  u  <  1/2  implies  that  the  leading  term  is  the  first, 

, -1  — 3  -1  | -3+u,  I  | 

i.e.  X  a  e  =  0(  x  )  as  x  -*•  °°. 


The  remaining  integral 


1-e 

l+e 

. 

_d_ 

I  dx 


iX (x  -1) 


p  (x,ax,t) 


ixe^ 

e  dx 


has  a  bound  of  the  form: 


_  , -4  .  5  -4  ^  4  -5  3  -6  2  -7  -8. 

CX  (a  e  +  a  c  +  a  e  +  a  g  +ae  ) 


which  has  leading  term 


[ixrrH 

as 

x  |  -*■<». 
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This  can  be  chosen  to  be  of  lower  order  than  the  contribution  from  the 


stationary  phase  points,  by  making 


7  9  5 

—  +  4u  <  -  —  ,  which  is  true  for  u  <  — 
2  4  16 


This  completes  the  proof  of  Lemma  3.5. 

To  finish  the  proof  of  Theorem  3.2,  we  remark  that  if  x  =  ct  +  6 , 


c  <  0  then  a  -*■ 


'  d  1/2 


12 


as  t  -*•  +°°,  and  the  usual  method  of  stationary 


-1/2 

phase  applies  [8] .  This  gives  a  decay  rate  of  X  which  is  proportional 


-1/2 

to  t  and  finishes  the  proof  of  Theorem  3.2. 


i 


1 


4.  GLOBAL  BEHAVIOR  AND  UNIQUENESS  FOR  THE  FIRST  ORDER  TERM  IN  THE 
N-SOLITON  WATER  WAVE  PROBLEM  . 

In  section  2  above,  a  formal  expansion  procedure  was  given  for  the 
Euler  equations  for  a  fluid  in  a  flat-bottomed  canal  which  was  near  the 
constant  horizontal  flow  of  Froude  number  1.  We  now  show  that  the  choice 
of  an  N-soliton  solution  of  the  KdV  equation  as  leading  term  in  this 
expansion  results  in  an  equation  for  the  first  order  term  which  has  a 
unique  "N-tuple  wave"  solution  if  we  add  a  symmetry  requirement.  As  noted 
previously,  this  term  satisfies  the  inhomogeneous  form  of  the  linearized 
KdV  equation: 


(4.1)  lu  ■  u,  +  u  -  6  (qu)  =  h(x,t) 

t  xxx  x 

where  cj(u,t)  is  an  N-soliton  and  h(x,t)  is  a  term  which  depends  only 
on  q(x,t).  A  simple  calculation  shows  that  in  fact  h(x,t)  is  a  linear 
combination  of  the  functions: 


a 


x  '  qqx  '  qxxx  '  q2qx  ' 


q  q  ,  q  q  ,  q 
^  ^xxx  X  XX  xxxxx 


We  remark  that  h(x,t)  contains  terms  of  the  form  F!(x,t),  which 

satisfy  the  linearized  equation.  It  is  rather  surprising  that  these 

"secular  terms"  [13]  do  not  give  rise  to  resonant  solutions.  The  usual 

choice  for  the  solution  to  Lu  =  F!  would  be  tF(  ,  which  grows  linearly 

ir,  t  in  the  moving  frame  in  which  =  x-4®,  t  remains  constant  as 

3  3 

t  *  However,  the  function  G_.' (x,t)  is  a  solution  of  the  homogeneous 
■  qua t ion  of  the  form  (see  the  discussion  in  the  Appendix) 


(4.2) 


G 1  ( x ,  t  )  =  c.  [  (x-12  62t)F’.  (x,t)  +  F  .  (x ,  t)  +  H1.  (x,t)  ] 


3 


3  3 


3 


SoS 


i  ■  a  rational  function  of  the  exponentials  (e  ^  which 
for  t  >  0  fixed.  Thus 
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(4.3) 


G ! (x,t) 

i-J  — — 


-  (x  43 jt)  F ! (x,t)  -  Fj (x,t) 


h!  (x,t) 


88  ^t  F ! (x,  t) 


and  since  L(G.!)=0,  we  have: 

(4.4)  L[(x-4Bjt)F*(x,t)  +  F^ (xf t)  +  H( (x,t) ]  =  L[86^t  F'(x,t)] 

=  8B*F»(x,t) 


The  function  (x-48..t)  f! (x,t)  +  Fj(x,t)  +  H((x,t)  has  the  property  that 
it  is  bounded  as  t  -*•  1=0  in  any  moving  frame,  even  £  .  =  constant,  so  we 
have  found  a  ' nonresonant'  solution  for  the  secular  term  F((x,t).  For 
the  secular  forcing  terms  G((x,t),  the  growth  in  the  obvious  solution  is 
quadratic  in  t  as  t  -*•  00  with  £ .  fixed  and,  to  the  best  of  our 
knowledge,  no  nonresonant  solutions  of  (4.1)  exist.  By  the  absence  of 

these  secular  terms,  the  perturbation  we  consider  is  rather  special. 

In  order  to  study  (4.1)  when  h(x,t)  is  a  linear  combination  of  the 

functions  listed  above,  we  use  the  following  representation  of  the  N-soliton 
solutions  of  the  KdV  equation  (see  Gardner,  Greene,  Kurskal,  and  Miura  [10], 
Tanaka  [22],  and  Deif t-Trubowitz  [7]): 


N 

(4.5)  q (x , t)  =  l  a.  F . (x,t) 

j=l  3  3 

2 

where  F^(x,t)  is  as  usual  the  squared  eigenfunction  f+ (x,  i  8.. ,  t)  .  Then, 
using  the  third  order  equation  satisfied  by  the  squared  eigenfunctions  [19], 
we  have: 


-35- 


(4.6) 


2 

q"'  =  l  a.[4(q(x/t)  +  B  )  F'  +  2q’F.; 
j=l  3  3  3  3 


=  6qqT  +  7  4a  .$  .F\ 

j=i  3  3  3 


2  N 

qq»  "=  6q  q1  +  £  4a .0 .q  F’ 

j=l  33  j 


q{5)  =  6qq»"  +  18 


q’q"  +  Y  [16a.B2(q+B2)F*.+8a.BZq’F. 

D  j  D  j  DID 


Thus  our  particular  forcing  term 
(4.7)  / F( (x, t) ,  q(x,t)F’ (x,t) , 


h(x,t)  is  in  the  span  of  the  functions: 
2 


q*  (x,t)F  (x,t) ,  q'q*  (x,t),  q»  q» * 


r  D  !».../ 


We  prove: 

Lemma  4,1.  Suppose  h(x,t)  is  a  linear  combination  of  the  functions  in 
(4.7).  Then  there  exists  a  solution  to  the  linearized  Kdv  equation  (4.1), 

Lu  =  h,  which  is  an  N-tuple  wave  in  the  following  sense: 

(i)  u(x,t)  -+  0  exponentially  fast  as  |x|  -*•  <*>  for  t  fixed 

(ii)  u(ct+<5,t)  ->■  0  exponentially  fast  as  t  -*•  +® 

2 

if  c  ^  46  ^  ,  j  =  1 , . . .  ,N 
2 

(iii)  lim  u  (4B.t+£rt)  exists  and  is  an  exponentially  decreasing 

tr«°  D 

function  of  £ 

(iv)  In  fact,  u(x,t)  is  a  sum  of  terms  which  are  either 

f  2 

(a)  rational  functions  in  the  N-exponentials  jexp  (Bj(x  -  4Bjt)} 

2 

with  the  same  denominator  [w(x,t)J  as  q(x,t)  _or 

(b)  of  the  form  (x-4B2t)  F((x,t). 

Proof  of  Lemma  4.1:  As  above,  we  write 


(4.8) 


L  (u) 


u  +  u  -  6  (qu) 
t  xxx  X 
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From  [10],  Theorem  3.6,  F((x,t)  satisfy  L(u)  =  0  and  F_.(x,t)  satisfy 

the  adjoint  equation 


(4.9)  v  +  v  -  6q  v  =0 

t  XXX  X 

Thus  L  (F  . )  =  -  6F.  q’  ;  L  ( tF '. )  =  F*  .  Also,  L  f  (x-12  g2t)  F '  +  F  .)  =  6q  FT  . 

3  3  3  3  v  3  3  3J  3 

2  2 

From  the  KdV  equation,  L(q  )  =  6q'q* f  -  6q  q1  and  L(q'')  =  ^q'q*'.  A  basis  fo 
the  solution  is  therefore  given  by 


(4.10) 

{ F  , 

tF],  (x-126  .  t) F ! +F  .  ,  q  ,  q1*} 

3 

3  3  3  3 

Using  L (G !) 

=  0 

and  (4.2),  an  equivalent  basis  is: 

(4.11) 

{v 

(x-462t)F»  H>,  q2,  q") 

l  3 

3  3  3  J 

The  functions 

Fy 

2 

q  ,  q11  have  properties  (i)-(iv) 

F  .  (x , t)  is  exponentially  decreasing.  In  the  Appendix,  we  show  that 
2 

(x-4Bjt)  F!  (x,t)  is  bounded  and  satisfies  (i)-(iv)  and  prove  that  H' (x,t) 
is  a  rational  function  of  the  expontentials  which  has  properties  (i)-(iv) 
as  well.  Assuming  these  results,  the  lemma  is  proved. 

Remarks:  (i)  Uniqueness:  If  we  choose  the  phases  of  the  N-soliton  so 

that  q(-x,-t)  =  q (x, t)  then  the  solution  in  Lemma  4  is  unique  provided 


we  require: 

(a)  u(-x,-t)  =  u(x,t) 

(b)  u(ct+6,t)  is  bounded  for  all  c  as  t  -*■  » 

2 

(c)  u(ct+6,t)  -v  0  exponentially  fast  if  c  ^  46 . . 

Proof :  From  the  results  of  Section  3,  the  kernel  of  L  is  spanned  by 

F^(x,t),  G*(x,t)  and  {  (  f 2)  '  (x ,k ,  t)  |  -oo  <  k  <  «>}  .  The  functions  F1. 

violate  (a);  G!(x,t)  violates  (fc)  for  c  =  46^;  by  stationary  phase 
2  2 

analysis,  (f  )'(x,k,t)  violates  (c)  for  c=-4k  (the  decay  is  algebraic, 
not  exponential).  Thus  u(x,t)  as  given  in  Lemma  4.1  is  unique,  since 
in  this  case,  the  functions  given  in  (4.11)  satisfy  all  these  conditions. 
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( 


1 

(ii)  Higher  order  terms:  Even  in  the  time-independent  case,  ! 

explicit  expressions  for  the  higher  order  terms  of  the  formal  expansion 
involve  more  complicated,  transcendental  functions.  For  even  the  second- 
order  term,  functions  like 

log  (1  - 

occur.  Thus  algebraic  methods  will 
like  T.,cmma  4.1  for  the  higher  order 


2 

tanh  8x)  sech  6x 

not  readily  yield  solvability  results 
terms. 


4 
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APPENDIX 


the  Wronskian  determinant  in  x  of  ip^, . . .  ,4/^. 

Then 

d2 

q(x,t)  -  -  2  — j  log  w(x,t)  . 
dx 

This  definition  is  sensible  because  w(x,t)  >  0  ,  which  we  show  below. 
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Lemma  A.l.  w(x,t)  >  0.  In  fact,  w(x,t)  is  a  sum  of  exponentials  with 


positive  coefficients. 


Proof:  w(x,t)  =  WN  cosh  ,  sinh  '  ***' 


SXp{BNCN}  +  (_1'  exp(”6NCN} 


W  ' 


By  the  multilinerarity  of  the  determinant,  this  is  the  sum  over  all  possible 


determinants , 


choices  e_.  =  ±  1  of  2  W^exp  e2  exp  ^2^2^ '  exp  ^3^3^' 

exp  ...  i.e.  we  have  upon  evaluating  the  Vandermonde 


(A. 2)  w(x,t)  = 


all  choices 

e  .  =  ±  1 
3 


{[l  c^cj}  e0cy 


2 [n/2] 


(£A  -  £jej)  • 


Since  0  <  <  *  *  *  <  '  the  num:ber  of  negative  factors  in  the 

product  can  be  found  explicitly.  Namely,  if  z.  , ...,e,  are  the  negative 

11  1r 

indices  for  a  given  choice  of  the  e's,  we  obtain  (i^-1)  +  •••  +  (i^-1) 

negative  factors  in  ]  T  (e,  8.  -  e.8.).  The  extra  factor  e„E  •••£__ 

j  £  k  k  j  j  24  2 [n/2] 

adds  an  additional  (-1)  factor  for  each  i^  which  is  even.  For  any 

choice  of  r  and  i^,...,i  ,  this  means  that  there  are  an  even  number 

of  (-1)  factors;  thus  every  term  in  the  sum  (A. 2)  has  a  positive 

N 

coefficient.  We  remark  that  all  exponents  £  z.&.£.  occur  in  w(x,t) 

j=l  333 

and  that  w(x,t)  =  w(-x,-t)  since  changing  {e,}  {-£_.}  does  not  alter 

the  coefficients  in  (A. 2).  This  proves  Lemma  A.l. 

The  eigenfunctions  f+(x,k,t)  are  given  explicitly  by: 


(A. 3)  f+(x,k,t) 


WN+1  AA'  ‘  ’  ’  A'  exP^ikx+4ik  t>) 
N  ~ ~ 

w(x,t)  ]  f  (ik  -  8  ) 

1=1 
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The  normalization  f  ^  exp(ikx+4ik Jtl  as  x  ■+  +00  for  t  fixed  is 

satisfied,  as  is  seen  by  looking  at  the  leading  term,  which  has  exponent 
N 

I  M-  (i.e.  pick  the  term  with  all  e  's  =  +  1  in  (A. 2)  and  the 
j=l  3  3  K 

corresponding  expansion  of  the  numerator).  From  the  fact  {[6],  [7])  that 
N 

T (k)  5  TT  (k+iB.)/(k-iB  .)  ,  we  obtain  a  similar  expression  for  f  (x,k,t) 
j=l  3  : 

using  T(k)  f_(x,k,t)  =  f+(x,k,t).  The  proof  that  f+(x,k,t)  satisfy  the 
Schrodinger  equation  with  potential  q(x,t)  defined  as  above  is  given  in 
Deift  [6];  the  basic  idea  is  to  use  Jacobi's  identity  for  the  Wronskians 
and  induction  on  N. 

From  (A. 3)  and  the  expression  for  f  (x,k,t),  it  is  easy  to  see  that 


f+(x,iB.,t)  +  (-1) j+1f_(x,iB .,t)  i  0  . 


Also,  from  (A. 2)  and  (A.  3),  it  is  clear  that  f+(x,ig_.,t)  decays  like 

exp  {-B.|£.|}  as  I Z,  .  j  00  since  the  exponentials  exp  {±8.5.1  in  the 
3  3  3  3  3 

numerator  will  cancel  each  other,  while  remaining  in  the  denominator  w(x,t) 
The  factor  is  g^(x,t)  defined  in  (3.6)  as 


(A. 4)  9 j  (x,t)  =  t-  (f_(x,k,t)  +  (-1)  3  +  1f+(x,k,t) ) 


k=i6. 


Differentiating  the  exponential  exp  { ikx+4ik  t}  gives  a  term 

2  o—i 

2 (x-128 jt) f+(x,iB j ,t)  ,  while  differentiating  the  factors  (±  ik) 

which  occur  in  the  (£,N+1)  entry  in  the  Wronskian  gives  terms  having 

N 

the  form  (c  exp{ X (x, t)) / (w (x, t))  where  X(x,t)  =  +6  .5  .  +  £  e  8.5,, 

^  3  J  £=1  ^  ^ 

whose  sum  we  denote  by  h^(x,t).  Thus  g^(x,t)  grows  like  exp  {8j|5j|l 
as  |5j|  00  and  is  of  the  form  2  (x-12P^t)  f+  (x ,  i8  ,  t)  +  h^(x,t)  where 

h^(x,t)  is  a  rational  function  of  the  N  exponentials  exp  {8^5^}  with 
denominator  w(x,t),  growing  like  exp  { 8 ^  |  ^  j  1 1  as  |  f;  |  ->  «>.  Since 
g^(x,t)  =  c^  f+(x,i6j,t)  g  j  (x,t)  ,  multiplying  by  f+(x,iB..,t)  we  have 
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Gj(x't)  2  2 

(A. 5)  - ~ — -  2 (x-12B jt) f+(x,iB j ,t)  +  f+<x,i8  j,t)h_. (x,t) 

e  2  (x-126 jt)  F.(x,t)  +  H.(x,t)  . 

Since  f^(x,if3j,t)  is  rational  with  denominator  w(x,t)  and  decays  like 
exp  t-3  j  (C  j  | }  as  _+  00 ,  H^(x,t)  is  rational  with  denominator 

(w(x,t))  and  is  bounded  as  00 •  Since  all  the  other  exponentials 

occur  in  the  numerator  with  growth  at  most  exp  {  £  2BJcJ}  and  these 

W 

terms  are  balanced  by  those  in  the  denominator,  H ^  (x,t)  is  bounded  for 
all  x,t  real.  It  then  follows  that  H.j(x,t)  is  a  sum  of  terms 

2 

which  decrease  exponentially  fast  as  t  “  except  in  the  frames  x-4B^t 

=  constant,  where  their  limit  is  an  exponentially  decreasing  function  of 

2 

the  variable  £  ■  =  x-4Bjt.  Note,  however,  that  for  t  -*■  00 ,  (x,t) 

"decouples"  into  N  exponentially  decreasing  bumps  moving  at  the  speeds 
4B£  with  the  same  phases  as  q(x,t)  as  t  ■+•  unlike  Fj(x,t),  these 
terms  give  rise  to  contributions  in  all  N  moving  frames.  These  are  the 
basic  properties  used  in  the  discussion  in  Sections  3  and  4  above. 
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